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I. INTRODUCTION 

Lepton asymmetry is not directly observable. Besides, 
in general, large lepton numbers L bigger than the baryon 
number B (L > B) and even L ^ B are not excluded by 
any profound theoretical principle. 

The origin of eventual big lepton numbers is not fully 
understood. There exist numerous different mechanisms 
for their generation. Large lepton numbers can be con- 
sistent with the observed small baryon number for ex- 
ample in the context of Grand unification, as pointed 
in [||; there exist different models producing L ^ B: 
in a scenario of baryogenesis with baryonic charge con- 
densate [||, where large and varying lepton asymmetry 
can be produced, the out-of equilibrium decays of heavy 
Majorana neutrinos can produce it as well, reso- 
nant neutrino oscillations can create considerable neu- 
trino asymmetry at low temperature See also the 
review papers . 

Besides, large lepton asymmetry has interesting im- 
li£ations in the early Universe: in baryogenesis models 
|l0[| , for solving the monopole problem ||ll[], as well 
as domain wall problem jl^ , in primordial nucleosynthe- 
sis jr^[T5|] for structure formation in the Universe [p^ , 
relic asymmetry can suppress neutrino transitions in the 
early Universe and have interesting implications for solar 
neutrino problem, as well as in cosmological nucleosyn- 
thesis (CN) 0. 

Therefore, in this work we consider it worthwhile to 
explore precisely the effect of lepton asymmetries and to 
study the problem of oscillations and their influence on 
primordial production of helium-4 in the presence of such 
asymmetries. Especially, in the present work we would 
like to show that much smaller asymmetries (|L| ^ 0.1) 
still can considerably affect cosmological nucleosynthesis, 
although indirectly through its effect on neutrino oscilla- 
tions. Thus, the presence of an initial asymmetry bigger 
than the baryon one may considerably change the bounds 
on oscillation parameters obtained from primordial nucle- 
osynthesis considerations. 

The theme of lepton asymmetry in nucleosynthesis 
models with oscillations was previously considered in sev- 



eral publications ||1 



The possibility of 



changing the nucleosynthesis constraints on lepton asym- 
metry due to oscillations was discussed within models 
with great L ^ 0.1 |2^. The dynamical evolution of L 
due to oscillations was studied for some particular values 
of L and neutrino mixing parameters |l8yi9| ]. The pos- 
sibility for loosening the nucleosynthesis bounds on mix- 
ing parameters due to the suppression of oscillations by 
great enough L was revealed . And also the possibility 
of generation of asymmetry due to resonant oscillations 
was studied 0,^, and some scenarios where oscillations 
generate asymmetry which further suppresses oscillations 
and as a result alleviates nucleosynthesis bounds on neu- 
trino mixing were discussed. However, always the sim- 
plifying assumption that the evolution of the neutrino 



ensemble follows the average neutrino momentum was 
made. And the case of high temperature resonant os- 
cillations Sm^ < with big mass differences and great 
asymmetries was considered. The nonresonant case was 
considered just for a certain set of parameters and for not 
very large initial L it was shown that neutrino asymmetry 
adjusts dynamically to zero ]l9t . 

In this work we consider the case of nonresonant oscil- 
lations Sm? > with small mass differences for a wide 
range of the initial values of L, namely 10~^° — 1, and 
provide a precise numerical study for the most unex- 
plored part of it: 10^*" < \L\ < 10~^. We insist here 
for the correct account for the spectral spread of neu- 
trino, which we find to play an important role in the 
oscillation-asymmetry interplay, and to be decisive for 
the realization of low temperature resonance in Sm? > 
case.| 

We find that even small asymmetries \L\ > 10~^ may 
considerably effect nucleosynthesis (contrary to the con- 
clusions of the previous works) and that asymmetry is 
also able to enhance oscillations and consequently to 
strengthen the nucleosynthesis bounds on neutrino mix- 
ing parameters for concrete regions in the parameter 
space {L,5m? ,{)), besides suppressing oscillations and 
thus loosening the bounds on mixing parameters from 
nucleosynthesis. The analysis of the asymmetry effect is 
provided using the precise kinetic approach to the prob- 
lem of neutrino oscillations and asymmetry, described 
in 0j2l[] . Actually, the present work is an extension of our 
previous study ||l|,0 for the case of a lepton asymmetry 
bigger than the baryon one. There we have calculated the 
bounds on the neutrino mixing parameters for a v,, ^ Vs 
mixing scheme with small mass differences Sm?' < 10"^ 
eV^, where is the active electron neutrino, while de- 
notes the sterile antineutrino. This case corresponds to 
nonequilibrium oscillations between electron neutrinos Ve 
and sterile neutrinos Vs when the latter do not thermalize 
till Ve decoupling at 2 MeV and oscillations become ef- 
fective after i/^ decoupling. We discussed the special case 
of nonequilibrium oscillations between weak interacting 
electron neutrinos and sterile neutrinos for small mass 
differences 5m^, as far as the case of large Sm^ was al- 
ready studied (both concerning small asymmetries L ^ B 
H-l^ and bigger ones L> B ||l|,|l^,§). Due to the cor- 
rect kinetic approach accounting for neutrino depletion, 
neutrino spectrum distortion and neutrino asymmetry for 
each neutrino momentum, we have obtained more precise 
bounds (an order of magnitude stronger than the existing 
previously) to the neutrino oscillation parameters. In the 
present work we will show how these bounds are changed 



*Mind that, as was shown in the pioneer work in models 
working with the mean momentum only, the resonance does 
not have place as far as actually the asymmetry calculated in 
this approximation reaches zero at resonance. 
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due to bigger lepton asymmetries. 

We would like to emphasize that the qualitatively new 
result, namely that small asymmetries may enhance oscil- 
lations and consequently to strengthen the nucleosynthe- 
sis bounds, was revealed only due to the correct approach 
accounting for the spectral spread of neutrino momenta 
and energy spectrum distortion due to oscillations, advo- 
cated in our previous papers |ll |2l] , |28|| and in the pioneer 
work of Dolgov ||2^ . Other studies missed this possibil- 
ity as far as there it was assumed that neutrino ensemble 
follows the behavior of average momentum even in the 
presence of oscillations and asymmetry. 

In this work we explore the different ways in which lep- 
ton asymmetry may affect oscillations (and vice versa - 
oscillations affect the asymmetry) and primordial nucle- 
osynthesis. We have defined numerically the concrete L 
ranges where these effects hold, namely: (a) big enough L 
- total suppression of oscillations corresponding to degen- 
erate nucleosynthesis, without oscillations (b) resonance 
region - L enhancing oscillations and leading to an over- 
production of He-4, and (c) small L - neghgible effect. 

The paper is organized as follows. In the following Sec- 
tion II we discuss neutrino degeneracy. In Section III we 
present the model of nonequilibrium neutrino oscillations 
and analyze the neutrino evolution, using kinetic equa- 
tions for the neutrino density matrix for each momentum 
mode. In Section IV we investigate <-> Vs the primor- 
dial production of helium in the presence of oscillations 
and neutrino degeneracy. The results and conclusions are 
presented in Section V. 

II. LEPTON ASYMMETRY - PRELIMINARIES 

At the epoch of interest - just prior to cosmological 
nucleosynthesis - at temperatures around few MeV, the 
lepton asymmetry is expressed through the asymmetries 
in the neutrino sector and the electron sector: 

L = Y,{N,-Ni)/N^, 

where £ = e, i^e, J^^i, ^^t, and i corresponds to their antipar- 
ticles. Big lepton asymmetry L > B may be contained 
only in the neutrino sector, as far as the electron-positron 
asymmetry must be equal to the proton-antiproton asym- 
metry from the charge conservation requirements. Hence, 
it will be more correct to talk about neutrino asymmetry 
or neutrino degeneracy, as it is most often referred to. 
It is convenient instead of {N^, — Nj^)/N^ to use the so 
called neutrino degeneracy parameter, ^ = ^/T, where /i 
is the corresponding neutrino chemical potential. In the 
thermal equilibrium of the early Universe at T = 2 MeV, 
the neutrino degeneracy for the case ^l/T <C 1 is given 
by 

i i 



where i = e, ^, t. 

There are no experimental or direct theoretical limita- 
tions for the magnitude or the sign of neutrino asymme- 
tries. However, there exist indirect limits on the neutrino 
degeneracies in the different neutrino flavors, obtained 
on the basis of astrophysical and cosmological considera- 
tions. Such are the constraints obtained from the present 
age and density of the Universe |Q ■ £, < 86 for the degen- 
eracy in only one neutrino species; from the requirement 
for long enough matter dominated period, necessary for 
successful structure formation: ^ < 6.9, and the strongest 
bound on the magnitude of L is obtained from primor- 
dial nucleosynthesis considerations ||l^,^,|l|: —0.06 < 
< 1.1 and l^^l < 6.9 for baryon to photon ratio 
7? < 1.9 X 10-^°. 

As far as primordial nucleosynthesis is particularly sen- 
sitive both to the neutrino degeneracies and to neutrino 
oscillations, it is interesting to consider the effect of neu- 
trino degeneracies in a modified model of primordial nu- 
cleosynthesis with neutrino oscillations. 

In this work we assume no degeneracy of muon and 
tau neutrinos for simplicity. However, this assumption is 
not essential. The results can be easily rescaled for the 
general case of asymmetry in all sectors. We consider 
the effect of neutrino asymmetries in the wide diapason 
of its initial values 10""'^° — 1. Both positive and negative 
values of L are considered. Asymmetries smaller than 
10^^*^ have negligible effect, therefore we have not stud- 
ied them. Now we have found that small asymmetries 
with magnitude \L\ < 10~^ have negligible effect^] also 
for the concrete mixing parameters of our model. We 
have proved that asymmetries greater than 10~* com- 
pletely suppress the oscillations and affect the nucleosyn- 
thesis due to the pure direct effects of neutrino degener- 
acy on nucleosynthesis - through affecting the expansion 
rate and the weak reaction rates. The latter effects have 
been studied already in the numerous publications on 
degeneracy and nucleosynthesis |^,^,^, therefore, we 
have not provided detail numerical calculations for such 
big asymmetries. While asymmetries in the intermediate 
region, namely 10~^ < \L\ < 10^* have not been studied 
systematically. 

In this work we prove that such asymmetries can both 
enhance and suppress oscillations depending on the con- 
crete set of mixing parameters. For that range we have 
provided detail numerical analysis of the effect of the neu- 
trino degeneracy on neutrino evolution and consequently 
on helium-4 production. Finally (as described in the fol- 
lowing sections) we have obtained the isohclium contours 
for different L values in the {Sm?,'d) plane and derived 
new limits on the neutrino mixing parameters for the 
cases with neutrino degeneracy. 



This confirmed our preliminary results stated in 
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III. EVOLUTION OF THE NEUTRINO 
ENSEMBLES IN THE PRESENCE OF INITIAL 
NEUTRINO DEGENERACY 

A. Neutrino oscillations - the model 

We suppose the existence of a sterile neutrino {SU (2)- 
singlet) i^s, and would like to explore the cosmological 
effect of nonresonant neutrino oscillations Ve ^ i^s on 
the primordial nucleosynthesis for a wide range of values 
of the initial lepton asymmetry. 

Within the model of interest oscillations proceed effec- 
tively after the active neutrino decoupling and till then 
the sterile neutrinos have not yet thermalizedQ so that 
their number density is negligible in comparison with the 
electron neutrino one. The model is described in detail 
in |,||. 

Oscillations between Vs (i^s = v^) and the active neu- 
trinos proceed according to the Majorana&Dirac [MhD) 
mixing scheme p2[ . For simplicity mixing present just in 
the electron sector is assumed Vi — Uu vi, I = e, s: 

Vi = CVe + SVs, 
V2 = -SVe + CVs, 

where Vg denotes the sterile electron antineutrino, c = 
cos(i5), s = sin(-!^) and ■& is the mixing angle in the elec- 
tron sector, the mass eigenstates vi and 1^2 are Majorana 
particles with masses correspondingly mi and m2. We 
consider the nonresonant case 5m? = m^ — fnl > 0, which 
corresponds in the small mixing angle limit to a sterile 
neutrino heavier than the active one. We would like to 
remind, that 'nonresonant' is used in the sense that it 
does not allow resonant transitions due to the finite tem- 
perature correction term (nonlocal term). However, as it 
is obvious, for the low temperature resonance this termi- 
nology is misleading. As far as for any sign of Sm^ then 
there is a possibility for a resonance due to the local term 
either in the neutrino sector or in the antineutrino sector. 



**This assumption of nonthermalization of steriles till elec- 
tron neutrino decoupling and the late effectiveness of oscilla- 
tions limits the allowed range of oscillation parameters for the 
discussed model: sin (2i9)(5m2 < 10"^ eV^ Q. However, it 
allowed us to examine the interesting nonequilibrium case of 
neutrino oscillations, in general with non Fermi-Dirac distri- 
bution, which is realized at such small mass differences. And 
besides, it allowed us to neglect the second order in Gf terms 
in the equations governing the neutrino evolution, and helps 
simplifying the calculation process. (For the nonequilibrium 
case studied, this is important, as far as accounting for the 
neutrino spectrum distortion for that case is obligatory, and 
it needs much longer calculation time.) 



B. The kinetics of nonequilibrium neutrino 
oscillations 

The kinetic equations for the density matrix of the 
nonequilibrium oscillating neutrinos in the primeval 
plasma of the Universe in the epoch previous to nucle- 
osynthesis, i.e. consisting of photons, neutrinos, elec- 
trons, nucleons, and the corresponding antiparticles, 
have the form: 



dgjt) 
dt 



Hp 



dgjt) 
dp 



+ ^[no,git)] + i[n^nt,g{t)], (1) 



where p is the momentum of electron neutrino and g is 
the density matrix of the massive Majorana neutrinos in 
momentum space. 

These equations account simultaneously for the partici- 
pation of neutrinos into expansion, oscillations and inter- 
actions with the medium. The first term in the equation 
describes the effect of expansion, the second is responsible 
for oscillations, the third accounts for forward neutrino 
scattering off the medium. The collisions (second order 
interactions in Gp of neutrinos with the medium) are 
neglected, as far as we consider the case of oscillations 
effective after neutrino decoupling. 

TCo is the free neutrino Hamiltonian in the eigenstate 
basis: 











while Tiint = a is the interaction Hamiltonian, where 

The first 'local' term in V is proportional to the 
fermion asymmetry of the plasma and for L > B is es- 
sentially expressed through the neutrino asymmetry: 



The second 'nonlocal' term is momentum dependent 
Q ^ T 1^;^ . When lepton asymmetries of the or- 
der of the baryon one are considered, the nonlocal term 
dominates at high temperature, while with cooling of 
the Universe in the process of expansion the local one 
becomes more important |2^Jl^,|l^. For greater asym- 
metries L > B the role of the nonlocal term decreases 
in comparison with the local one, and usually for tem- 
peratures less than a few MeV it is neglected. We have 
precisely accounted for it in this work, having in mind 
the following. Due to the momentum dependence of the 
nonlocal term even when for the mean neutrino momen- 
tum modes at a given temperature this term is smaller 
than the local one, for the high energetic modes it may 
still play some role. We have checked numerically that 
the nonlocal term can change from a few up to about f 4% 
the results on helium-4 production. The maximum effect 
being, as can be easily guessed, in the case of maximum 
mixing. This is noticeable. 
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An analogues equations hold for the antineutrino den- 
sity matrix, the only difference being in the sign of the 
fermion asymmetry: £ is replaced by — £. 

Medium terms depend on neutrino energy density 
and neutrino degeneracy, thus introducing a nonlin- 
ear feedback mechanism. As was already observed in 
previous works oscillations change neutrino- 

antineutrino asymmetry and it in turn affects oscillations. 
For large L the evolution of neutrino and antineutrino 
ensembles is strongly coupled and hence, it must be con- 
sidered simultaneously. 

The equation results into a set of coupled nonlinear 
integro-differential equations with time dependent coef- 
ficients for the components of the density matrix of neu- 
trino. Due to the greater L the coupling is stronger than 
for the case L = 10~^°. As far as for these strongly cou- 
pled nonlinear equations an analytic solution is hardly 
possible, we have provided an exact kinetic analysis of 
the neutrino evolution by a numerical integration of the 
kinetic equations for the neutrino density matrix for each 
momentum mode. The evolution of the neutrino ensem- 
bles is followed numerically from the freezing at 2 MeV 
till the formation of helium-4. 

Our numerical analysis shows that the coupling be- 
tween the neutrino and antineutrino ensembles leads to 
their similar behavior: Whenever the resonance condi- 
tion is fulfilled for neutrino (or antineutrino) , and the en- 
semble suffers a resonant oscillations, due to the strong 
coupling between the systems, the antineutrino ensem- 
ble shows the same behavior (after some negligible delay 
time) too.[^ The results have almost negligible depen- 
dence on the sign of the initial asymmetry, as could be 
expected from this behavior of the ensembles. The con- 
crete 'spectral' mechanism of that strong coupling be- 
tween the ensembles in our case is described in detail 
further on. 

The initial condition for the neutrino ensembles in the 
interaction basis is assumed of the form: 



1 




where = exp(e - E^/T)/{1 + exp(e - E^/T)). 

It corresponds to an equilibrium distribution of active 
electron neutrinos, and an absence of the sterile one. We 
have analyzed the evolution of nonequilibrium oscillat- 
ing neutrinos by numerically integrating the kinetic equa- 
tions for the period after the electron neutrino decoupling 



Opposite to the traditional naive conclusions, based on 
observations just for a certain parameter values and provided 
for the mean neutrino momentum, that in the low tempera- 
ture resonance case (when the nonlocal term is neglected in 
comparison with the local one) the resonance may take place 
either in the neutrino sector or in the antineutrino one, but 
not both [Fi. 



till the freeze out of the neutron-proton ratio (n/p-ratio), 
i.e. for the temperature interval from 2 MeV till 0.3 MeV. 
We have explored the problem using the Simpson method 
for integration and the fourth order Runge-Kutta algo- 
rithm for the solution of the differential equations. In 
order to keep the appropriate precision of calculation for 
the wide range of initial parameter values, we used a step 
for the evolution with temperature decrease of the order 
10^^ — 10~^ MeV and for the calculation of the neutrino 
spectrum a step for E^/T of the order 10~^ — 10~^. 

The oscillation parameters range studied is 5vn? G 
[10"", 10"^] eV^ and -d £ [0, 7r/4], which is the full set of 
mixing parameters in the discussed model |^ of oscil- 
lations. The neutrino degeneracy was varied for a wide 
range 10^^° — 10"'' of L^, as discussed in the previous 
chapter. 

The distributions of electrons and positrons were taken 
the equilibrium ones. The neutron and proton number 
densities, used in the kinetic equations for neutrinos, were 
substituted from the numerical calculations in CN code 
accounting for neutrino oscillations. I.e. we have simul- 
taneously solved the equations governing the evolution of 
neutrino ensembles and those describing the evolution of 
the nucleons (see the next section) . The baryon asymme- 
try /3, parameterized as the ratio of the baryon number 
density to the photon number density, was taken to be 
3 X 10-1°. 



IV. NUCLEOSYNTHESIS WITH 
NONEQUILIBRIUM OSCILLATING NEUTRINOS 

We have calculated the production of helium-4 in a 
detail model of primordial nucleosynthesis, accounting 
for the direct kinetic effects of oscillations and neutrino 
degeneracy on the neutron-to-proton transitions. The 
effect of oscillations on helium-4 has been discussed in 
numerous publications 34- |3^ , p"^ -^. We follow the ap- 
proach described in ||l|,21|. The numerical analysis is 
performed working with exact kinetic equations for the 
nucleon number densities and neutrino density matrix 
in momentum space. This enabled us to analyze the 
direct influence of oscillations onto the kinetics of the 
neutron-to-proton transfers and to account precisely for 
the neutrino population depletion, distortion of the neu- 
trino spectrum and the role and dynamical evolution of 
neutrino-antineutrino asymmetry. 

The master equation, describing the evolution of the 
neutron number density in momentum space n„ for the 
case of oscillating neutrinos I'g <-> z/s, reads: 

(dUn/dt) = Hpn {dUn/dpn) + 

+ j df7(e",p, v)\A{e^ p vn)\'^ 
dn{e+ ,p,C')\A{e+n pD)\^ 
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X [ne+n„(l - Qee) " TlpQeei^ " "e+)] 



(2) 



We have calculated the neutron-to-photon ratio evolu- 
tion prior to nucleosynthesis epoch, namely for the tem- 
perature interval T E [0.3, 2.0] MeV and for the full range 
of oscillation parameters of our model and L in the range 
10-1° - IQ-''. The value of Qee{Eu/T) at each integra- 
tion step was taken from the simultaneously performed 
integration of the set of equations for neutrinos, i.e. the 
evolution of neutrino and the nucleons was followed self 
consistently. 

The initial values dXT — 2 MeV for the neutron, pro- 
ton and electron number densities are their equilibrium 
values. The parameters values of the CN model, adopted 
in our calculations, are the following: the mean neutron 
lifetime is r = 887 sec, which corresponds to the present 
weighted average value, the effective number of relativis- 
tic flavor types of neutrinos during the nucleosynthesis 
epoch is assumed equal to the standard value 3. 

On the basis of this analysis the primordially produced 
He-4 value was obtained as a function of neutrino degen- 
eracy parameter and neutrino mixing parameters. 



V. RESULTS AND CONCLUSIONS 

We have calculated the evolution of the neutron-to- 
nucleon number densities Xn{t) — Nn{t)/{Np + Nn) for 
each set of values in (L,6m'^,'d) parameter space. In 
Figs. 1 the evolution of X„ is illustrated for initial asym- 
metry of the order of the baryon one and for L — 10"^ 
and L = 10"^. As is seen from the figures, even relatively 
small asymmetry (L = 10"^) can considerably effect the 
neutron-to nucleon ratio. In some cases, i.e. for maximal 
mixing (Fig. la) leading to underproduction of helium-4, 
while in others (Fig. lb) to its overproduction in compar- 
ison with helium yields calculated in models with small 
asymmetry [L ^ 10" i°). 

In Fig. 2 the dependence of the frozen neutron num- 
ber density relative to nucleons X^{L) on the the initial 
neutrino asymmetry for sin^(2i?) = ^^^^ different 

fixed Sw? is illustrated. The dependence of the frozen 
neutron number density relative to nucleons X^{'d) on 
the mixing angle for different fixed (5m^, is presented in 
Fig. 3. 

It is interesting to note, that, as illustrated in Fig. 3 
depending on the mixing angle, it is possible to weaken 
or to strengthen the limits in comparison with the model 
of nucleosynthesis with oscillations without a consider- 
able lepton asymmetry. Namely, for large mixing, the 
presence of neutrino degeneracy of the order 10"^, for 
example, leads to a less strong constraints on the mixing 
parameters, weakening the oscillation effect of overpro- 
ducing helium. While for smaller mixing angles the same 
asymmetry leads to stronger limits on mixing parame- 
ters, as a result of the enhancement of oscillations, and 
hence increasing the helium overproduction. 



These results can be understood as follows: As it was 
shown in |^l| in case of baryon like small lepton asym- 
metry the kinetic effects (neutrino population depletion 
and distortion of neutrino spectrum) due to oscillations 
play an important role and lead to a considerable over- 
production of helium. 

In the presence of a bigger asymmetry both the neu- 
trino depletion and spectral distortion are changed due 
to the interconnections between oscillations and asymme- 
try. Large enough asymmetries (in our case L > 10"^) 
affect oscillations, while strong enough oscillations also 
influence asymmetry and lead to its dynamical evolu- 
tion. We would like to note, that even for relatively 
small \L\ < 10^^ the effect of asymmetry is consider- 
able (contrary to the thought before). And also that 
even oscillations with small mass differences are able to 
influence neutrino asymmetry. A characteristic situation 
in our model is that great asymmetries cannot be cre- 
ated: neutrino asymmetries typically of an order of mag- 
nitude greater than the initial one are generated, however 
they are rapidly oscillating, the average value being zero. 
In our case of nonresonant oscillations, the presence of 
the second system (for antineutrino) hinders the contin- 
uous growth of asymmetry. In more detail asymmetry- 
oscillation relations will be analyzed in a following paper. 
Our main purpose here was the study of asymmetry ef- 
fect on primordial production of helium-4, and derivation 
of the nucleosynthesis limits on neutrino mixing in the 
presence of neutrino degeneracy. 

For the oscillation model discussed we have found three 
interesting regions for the range of the lepton asymmetry. 
To be more clear let us consider nearly maximal mixing: 
sin2(2i5) = 10-0-05 (see Fig. 2 for illustration). Then 
these regions are as follows. 

A. In the range \L,j\ > 5 x 10^^ the asymmetry fully 
suppresses neutrino oscillations .|^ The yield of helium-4 
coincides with the values obtained in the models of pri- 
mordial nucleosynthesis without oscillations. Therefore, 
the nucleosynthesis bounds on the mixing parameters in 
the presence of such large asymmetries are waved away. 

For \L\ < 10"^ the direct kinetic effect of such asymme- 
tries on the neutron-to-proton transfer is almost negligi- 
ble, i.e. for the range 10"'' — 10"^ the only role of asym- 
metry is to suppress oscillations. Asymmetries greater 
than that effect nucleosynthesis by directly changing the 
kinetics of nucleons and are exhaustively studied in pre- 
vious works ||l^,|lj,0 . And the effect of the asymmetry 
is well known: Neutrino degeneracy effects element pro- 
duction in two ways |13|. First, the nonzero chemical 
potential results to an increase of neutrino energy den- 
sity and hence to an increase of the expansion rate of 



"'"'"This result is in accordance with the estimations made in 
[ p^ , where it was shown that a big neutrino asymmetry may 
lead to suppression of oscillations. 
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the Universe thus allowing less time for the nuclear re- 
actions to proceed. Second, it alters the reaction rates 
governing the neutron-proton transitions, as far as these 
rates are extremely sensitive to the neutrino spectrum. 
Namely, electron neutrino degeneracy results into pre- 
ponderance of protons, hence into a drop of primordial 
helium-4 abundance, while electron antincutrino degen- 
eracy leads to the increase of neutrons. Therefore for 
L > 10~^ we have not provided detail numerical calcula- 
tions. 

B. The range of smaller than 10^^ asymmetries was 
more appealing to us as far as it was totally unexplored 
till now. Therefore, we numerically analyzed the problem 
for the initial values of the neutrino asymmetry in the 
range 10"^" — 10~^, and for the full mixing parameter 
space of the active-sterile oscillations model, described 
in i.e. for all mass differences Srn^ < 10~^ eV^ and 
mixing angles i!). 

The zone 10"'^ < |L| < 5 X 10"^ is the most interest- 
ing one. The asymmetry with magnitudes of that order is 
strong enough to influence neutrino oscillations. And al- 
though the neutrino asymmetry is too small to effect nu- 
cleosynthesis directly, it does effect nucleosynthesis yields 
in indirect way through its influence on neutrino oscilla- 
tions pattern. The asymmetry values in that range, de- 
pending on the concrete values of mixing parameters, are 
able to enhance oscillations. So, depending on the mixing 
angle for some L in this interval the ensemble of neutri- 
nos expires a resonance. For nearly maximal mixing it 
is roughly around L ~ 10^^ (Fig- 2). This enhancement 
of oscillations is big enough to influence considerably the 
electron neutrino and electron antineutrino number den- 
sity and spectrum i.e. it leads to an enhanced depletion 
of the number densities of neutrinos and antineutrinos(!) 
and a decrease in their mean energy, as well as nontrivial 
evolution of the asymmetry itself. Hence influencing the 
helium-4 primordial value, which is extremely sensitive 
to it. The result is overproduction of helium-4 in com- 
parison with the case of primordial nucleosynthesis in the 
presence of oscillations with negligibly small L. This on 
its turn leads to a strengthening of the nucleosynthesis 
bounds on the neutrino mixing parameters for that cer- 
tain set of parameters (L, Sm"^, d). 

The ability of L within that range to enhance the oscil- 
lations looks quite amazing at a first glance. The naive 
picture one expects as a result of increasing the initial 
L is a gradual suppression of oscillations proportional to 
L. As far as L value, calculated for the mean neutrino 
momentum is by orders of magnitude bigger than the 
necessary one for a resonance transfer L ^ Lr. 

However, the detail numerical analysis, accounting for 
the momentum spread of neutrino ensemble, showed a 
more complex picture. Varying L for fixed mixing pa- 
rameters we have observed a resonance region, i.e. an 
enhancement of oscillations, as seen in Fig. 2 and Fig. 3. 
And besides that, when a resonance for a given L value 
was observed for neutrino it was followed with a negli- 
gible time delay by a resonance in antineutrino! This 



leads to a considerable depletion of number densities and 
to further overproduction of helium-4. This amazing re- 
sults can be explained as follows. 

The system of nonlinear differential equations cannot 
be solved analytically without radical assumptions. How- 
ever, the qualitative behavior of the ensembles and the 
obtained results concerning helium-4 can be guessed from 
the following simplified considerations. 

The resonant conditions for the neutrino looks like : 

cos{2^){dm'^/2E^) = V2Gf{L - Q/M^)N^ 

while for antineutrino the sign before L term is the op- 
posite one. 

For the mean neutrino momenta p = 3.15 T in our 
model 

|L|» Q/M^ 

|L|> coB{2d)5m'^ I{2^/2E^GfN^). 

Now let us consider the possibility that at a given tem- 
perature, for neutrino of a given momentum p < p the 
resonance condition is fulfilled, i.e. L has the resonance 
value Lr{p). Then neutrinos with this momentum suffer 
a resonant transfer, leading to a decrease in the number 
densities of neutrinos (in favour of the sterile neutrinos). 
As far as for antineutrinos L has the opposite sign, oscil- 
lations remain suppressed, and the number densities of 
antineutrinos do not change. Hence, the net result is a de- 
crease in L due to this resonant transition. This decrease 
makes possible the fulfillment of the resonant condition 
for more energetic neutrinos, leading to further decrease 
of L and so on; due to this 'resonance wave' passing to 
neutrinos with higher and higher momenta, the neutrino 
number densities expire a considerable depletion, conse- 
quently, having in mind the small initial values of the 
neutrino asymmetry {SN^, ~ L^Ny) soon this resonant 
wave leads to a change in the sign of L. This suppresses 
further resonant transfer for neutrinos, however, for the 
antineutrino ensemble now there appears the possibil- 
ity for a resonance. So, the same process follows in the 
antineutrino ensemble. The only difference being that 
contrary to the neutrino system, for antineutrinos the 
resonance is first fulfilled for the high energy neutrinos 
and then passes to the low energetic ones - i.e. it reaches 
more rapidly neutrinos with mean momentum, and the 
process is more avalanche like. This rapid decrease of 
antineutrinos leads to a rapid bump of L which again 
becomes positive. This pendulum like process proceeds 
relatively fast, the oscillation of the asymmetry proceeds 
much faster than the temperature decrease. The total 
effect is that the resonant transfer both for neutrino and 
antineutrino system is realized even for a considerable 
initial asymmetry values (nonresonant ones, when calcu- 
lated for the mean neutrino momentum), and due to the 
enhanced transfer of active to sterile neutrinos a consid- 
erable 'resonant' production of helium-4 is observed. 

C. For small asymmetries \L\ < 10~^ we have con- 
firmed our preliminary results from [pl|, namely that 
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asymmetry with these values has too week effect on nu- 
cleosynthesis to be considered. 

Finally, we have calculated the He-4 dependence on 
the neutrino asymmetry for the whole range of mixing 
parameters of our model and wc have obtained the pri- 
mordial helium yields Yp{L,Sm^,i!}). This enabled us to 
obtain isohelium contours in the Srn^ — d plane. Some 
of these constant helium contours, for L — 10"^° and 
L = 10^^, are presented in Fig. 4. 

On the basis of these results, requiring an agreement 
between the theoretically predicted and the inferred from 
observations values of primordial helium, we have ob- 
tained cosmological constraints on the neutrino mix- 
ing parameters corresponding to different initial lepton 
asymmetry values. Assuming the conventional observa- 
tional bound on primordial ^iJe Yp = 0.24 the cosmo- 
logically excluded region for the oscillation parameters 
on the plane sin^(2??) — Sm^ in Fig. 4 lies to the right 
of the Yp = 0.245 curve, which gives 5% overproduction 
of helium in comparison with the standard value. For 
comparison the curves corresponding to i = 10"^" and 
L = 10"^ are plotted. 

In conclusion we would like to stress once again that 
lepton asymmetries \L\ > 10~^ can both enhance or sup- 
press (completely or partially, depending on the concrete 
values of the model parameters) neutrino oscillations and 
consequently can strengthen or weaken oscillation effect 
on primordial nucleosynthesis. Therefore, the constraints 
on the oscillation parameters obtained from nucleosyn- 
thesis with lepton asymmetry of the order of the baryon 
one ]l8|-p6|,p| can be considerably changed - either relaxed 
(or even totally removed) or strengthened for (note!) not 
very large values of the initial lepton asymmetry. This 
may have interesting astrophysical and cosmological im- 
plications. 

Vice versa, having definite values for the neutrino mix- 
ing parameters, it will be possible, on the basis of require- 
ment of agreement between the theoretically calculated 
and the extracted from observations values of helium-4, 
to put cosmological constraints on L even within that 
range of exclusively small magnitudes. 

Finally we would like to point to the essential dif- 
ferences between our work and the previous ones: we 
have considered the nonresonant case of oscillations with 
very small mass differences; the oscillations proceed ef- 
fectively relatively late i.e. after neutrino decoupling and 
the change in the lepton number is dominated by oscilla- 
tions, not by collisions; wc have accounted precisely for 
the momentum spread of neutrinos. 

Having in mind the latest news from SuperKamioka 
concerning the solar neutrino oscillation solutions ^T\ , 
the very small mass differences concerned in our paper 
may look much more attractive, and hence this study 
will be interesting for the solar neutrino audience too. 
On the other hand we would like to stress that mod- 
els exploiting the suppression of oscillations due to L or 
the generation of L which further suppresses oscillations 
must be exploited with caution and be carefully calcu- 



lated for the concrete {L,5vn?,d) values of interest. As 
far as it was shown here that the effect of asymmetry is 
not straightforward. Asymmetry - oscillations interplay 
reveals greater complexity than the rough estimations 
suggest and should be studied carefully. 
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Figure Captions 



Figure la: The evolution of the neutron number density relative to nucleons = Nn{t)/{Np + Nn) for the 

case of oscillations with maximal mixing and Sm^ = 10~^ eV^ for different values of the initial asymmetry { L = 10~^, 
L = 10-5 and L = 10-i°) is plotted. 

Figure lb: The evolution of the neutron number density relative to nucleons Xn{t) = Nn{t)/{Np + Nn) for the case 
of oscillations with sin^(2'!?) = iQ-^ o^ and Sm^ = eV^ for different values of the initial asymmetry { L = 10~^, 
L = 10-5 and L = IQ-iO) is plotted. 

Figure 2: The figure illustrates the dependence of the frozen neutron number density relative to nucleons X„ = 
Nn/iNp + Nn) on the value of the initial neutrino asymmetry for sin^(2-i?) = lO-^'''^, and different mass differences. 
For comparison the standard curve is plotted also with dashed line. 

Figure 3: The dependence of the neutron to nucleon freezing ratio on the mixing angle for L ~ IQ-^ for different 
mass differences 8rr? = IQ- eV^ and 6vp? = 10-* eV^ is shown. For comparison with dashed lines the corresponding 
curves with small asymmetry L = 10-^° are presented. 

Figure 4: On the dm?' — d plane the constant helium contours calculated in the discussed model of cosmological 
nucleosynthesis with neutrino oscillations for L = 10-^ and L = 10~^° are shown. 
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Figure la: The evolution of tlie neutron number density relative to nucleons = Nn{t)/{Np + Nn) for the 

case of oscillations with maximal mixing and Sm^ = 10~^ eV^ for different values of the initial asymmetry { L = 10~^, 
L = 10-5 and L = 10-^°) is plotted. 




Figure lb: The evolution of the neutron number density relative to nucleons Xn{t) = Nn{t)/{Np + Nn) for the case 
of oscillations with sin^(2i9) = lO"^ ^^ and Sm'^ = 10~^ eV^ for different values of the initial asymmetry { L = 10~^, 
L = 10-5 and L = 10-^°) is plotted. 



log(sin 2i^) = -0.05 

0.144 n 




Figure 2: The figure illustrates the dependence of the frozen neutron number density relative to nucleons X„ = 
Nnl{Np + Nn) on the value of the initial neutrino asymmetry for sin^(2i?) = 10"° '^^, and different mass differences. 
For comparison the standard curve is plotted also with dashed line. 




Figure 3: The dependence of the neutron to nucleon freezing ratio on the mixing angle for L = 10~^ for different 
mass differences Sm^ = 10~^ eV^ and Sm^ = 10~* eV^ is shown. For comparison with dashed lines the corresponding 
curves with small asymmetry L = 10~^° are presented. 




Figure 4: On the Sni^ — i? plane the constant helium contours calculated in the discussed model of cosmological 
nucleosynthesis with neutrino oscillations for L = 10~^ and L = 10~^° are shown. 



